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We study some examples of collective behavior in neutrino clouds governed by the neutral-current 
neutrino-neutrino interaction. The standard equations for analyzing such systems are rederived 
in a two-step process: first, a replacement of the full interaction Hamiltonian with a "forward" 
Hamiltonian that contains only the momentum states that were initially occupied by a neutrino of 
one flavor or another; second, a factorization assumption that reduces the time evolution problem 
to the solution of coupled equations for the expectations of various bilinear forms in the neutrino 
fields. We designate the latter as the "classical" equations. We analyze some solutions of these 
equations in cases in which the initial momentum and flavor distributions of neutrinos are strongly 
anisotropic in space. In some cases we find an instability that leads to rapid evolution of the flavor- 
angle distribution, even when it is seeded by a very small initial flavor mixing (or alternatively 
by a very small neutrino mass 2 term). Turning to the more complete case in which we do not 
assume the classical factorization, but instead solve for the evolution under the influence of the full 
"forward" Hamiltonian, we find the possibility of rapid evolution, under our definition, even when 
there is no seeding from conventional neutrino mixing. This "speed-up", which occurs in exactly 
the same parameter range as do the instabilities in the classical case, to some degree confirms earlier 
conjectures of such collective speed-up's. However the time scales found in the present work are 
larger than those previously conjectured time scales by a factor of order logiV, where N is the 
number of neutrinos in a volume of dimension of the reaction length. Even with this lengthening of 
the "speeded-up" scale there can be situations in which these effects will dominate the short term 
behavior of a system. 
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1. INTRODUCTION 

Clouds of neutrinos, within media that may or may not 
have other constituents, are capable of behaving in many 
different interesting ways, some of which may even be 
relevent to real astrophysical situations. Within the con- 
text of the theory in which the neutrinos have standard- 
model couplings to baryons and leptons, augmented by 
a flavor-mixing mass term, we enumerate some possible 
behaviors that go beyond vacuum neutrino oscillations: 

1. The familiar index of refraction or "matter" effect 
from forward scattering on electrons 0], which 
enters the interpretation of the results of some solar 
neutrino experiments. 

2. Slowing of oscillation rates in the presence of 
"flavor-measuring" noise Q-jl]. 

3. Synchronization of oscillations in the presence of 
"flavor-blind" noise, [6j. 

4. Synchronization of oscillations through the nonlin- 
ear effects of the neutral-current neutrino-neutrino 
interaction 0-0- 

■5. Many-body effects that may 15], 16], or may not 
[l7|. speed-up some neutrino reaction rates. 
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In the present paper we discuss yet another set of be- 
haviors, closely linked to items #4 and #5, above, but 
with potential physical effects that are distinct. To set 
the context we note three different time-scales that can 
enter a system in which only neutrinos are present. 

a) Iscat = (G"pE 2 n v )~ 1 , the characteristic time for 
ordinary scattering of a neutrino in a cloud, where E is 
the energy scale of the neutrinos in the cloud and n„ is 
the neutrino number density. In what follows we consider 
only times that are much less than T scat . 

b) Tg = (Gpn^ 1 , the time over which the phase of 
the wave function of a neutrino is changed by an amount 
of order 2?r, due to interactions with other neutrinos in 
the cloud. We refer to T$ as the "short time scale" . 

c) Tl = (dm 2 /E)^ 1 , the neutrino oscillation time 
scale, which we designate as the "long time scale" . 

Under all conditions that we shall consider, which in- 
clude the parameter domains for both early-universe and 
supernova applications, we will have T scat >> Tl >> Tg. 
In this case, the evolution over short time scales will be 
determined by a "forward" neutral current interaction 
Hamiltonian, which we define in detail later. The evo- 
lution over the neutrino oscillation time-scale is deter- 
mined by a combination of the neutrino-mass terms and 
the neutral-current terms. These two effects are com- 
bined in a nonlinear equation for the neutrino distribu- 
tion function that has been presented in refs. 0], 0] 
and elsewhere and has been discussed and applied by nu- 
merous authors 201- 241. In certain circumstances the re- 
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suit is the synchronization of neutrino oscillation modes, 
as mentioned in # 4 in the above list. In the appli- 
cations of these nonlinear evolution equations, whether 
or not synchronization is operative, authors have usually 
assumed an isotropic medium. Because the forward scat- 
tering of neutrinos is strongly dependent on the angle of 
incidence, the equations for the anisotropic case are more 
complex. 

Here we study some examples of violently non-isotropic 
cases, and find behavior that is qualitatively very dif- 
ferent from that in isotropic cases. For example, in 
the isotropic case, if the system begins in a near flavor- 
eigenstate, with a small flavor mixing as a perturbation, 
then there will ensue slow oscillatory behavior, with a 
time scale Tj,. In some non-isotropic examples, on the 
other hand, we will find exponential growth of the per- 
turbation, on the short time scale, T$, followed by very 
non-linear and large amplitude oscillations, still on the 
short time scale. These oscillations do not change the 
total numbers of neutrinos of a particular flavor in a sys- 
tem, nor the all-over energy distribution or angular dis- 
tribution, but they can quickly change the angular and 
spectral distributions of a particular flavor. When we ex- 
plore, to some degree, the space of combined anisotropy 
and initial arrangement of flavors, we find that there are 
stable and unstable regions and that the boundaries will 
be rather hard to classify. These effects are surely of in- 
terest with respect to supernova physics, and could con- 
ceivably be of interest with respect to deviant theories of 
early universe evolution. 

Following many previous studies, we begin with a sys- 
tem that contains two flavors of neutrinos and with equa- 
tions that capture the complete effects of forward scat- 
tering. We make our considerations absolutely concrete 
by beginning with a "totally forward" effective Hamilto- 
nian, iff or , which contains both neutrino oscillations and 
the neutral current interactions. We expand the neutral- 
current Hamiltonian in terms of the creation and annihi- 
lation operators (a£)t, <x£, (af af for neutrinos of flavors 
r and s. Taking the set of momentum states {pi} that 
are occupied in the initial state, by neutrinos of cither 
species v r , v s (or mixtures thereof) and dropping from 
H all of the operators referring to other momenta, we 
obtain 

H ioT = ■ 3i + — =— Y <k,jai ■ <jj + H R , (1) 

i i,j 

where the operators Si are given by, 

= (atfa* + (aS)ta; , 
of =-i{a^a s +i{a^a] , 

of = (ajjtaj - ( a «)t $ . (2) 

Using the usual anticommutation rules for the operators 
a, yields the commutation rules of Pauli spin operators 



for the <7j. The coefficients djj are given by, 

(k,j = 1 -Pi -Pj/(\Pi\ Wj\). (3) 

The neutrino oscillation parameters are contained in A,; = 
Sm 2 1 pi, providing the scale, and in the vector B, which 
we take to have unit length. 

The residual part of QJ, Hr, contains only terms that 
commute with all of the operators a!- , these consist 
of the kinetic energy terms and all terms involving only 
the operators (ap^aj + (a|)ta|; the contribution of these 
terms to the amplitudes that we calculate will be all-over 
phase factors. We will occasionally use the language of 
spins on sites to characterize the states of the system 
governed by QJ, while understanding that the sites are 
really momentum states and the spins are really flavors. 
We write the Heisenberg equations of motion for the op- 
erators Si(t), 

-^■CTj = 2\iB x di - ^£jL ^2 di,j Si x & j . (4) 

3 

All of the results reported in refs. ^3"H3 depend on 
the following two approximations : 

(a.) Replacement of the equation of motion Q by its 
expectation value, (), with the assumption that, on the 
RHS, we can make the replacement, 

{ffi x B 3 ) -> {Si) x (Sj) , (5) 

as though the spins on different sites are not correlated. 
We shall refer to the resulting equations as "classical" ; 
having used the Heisenberg commutator to derive the 
equations, we then treat the expectation values as clas- 
sical variables. In what follows, we take the initial wave 
function of the system to embody the assumption J5j. 
But does not continue to hold as time progresses. 
Arguments have been made in the literature, albeit im- 
plicitly, that in the limit in which the number of parti- 
cles N — > oo equations like JSJ will be applicable over 
a large time scale, classical behavior being expected for 
large "spins" . We shall examine the validity of this as- 
sumption later. 

(b.) Isotropy of all distributions. This appears to allow 
us to make the replacement djj — > 1. At the very least, 
the applicability of this assumption depends on the ap- 
plication that is being pursued. Clearly, and as noted, e. 
g., in ref. j^, supernova applications require the reten- 
tion of the angular dependence embodied in the complete 
dij 's. In (nearly) isotropic early-universe applications it 
is probably adequate to use d^j — > 1, as in refs. [2l| : 
however, it is conceivable that instabilities of non-istropic 
distributions in the (nearly) collisionless neutrino sector 
could combine with the fluid dynamics of the collisional 
component of the matter to produce interesting phenom- 
ena. 
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In the main body of this paper we shall adopt assump- 
tion (a.), represented by (J5J, but drop the assumption 
(b.) of isotropy. We return in sec 6 to examination of 
the limitations of assumption (a.). 

In fllf. the operators, cr^, act on the wave functions 
of individual particles and there is an extensive factor of 
volume. To introduce operators with equations of motion 
and initial conditions that involve intensive variables, or 
densities in the present case, we divide the solid angle for 
momentum directions into a set of regions {a, b...} each 
with a reasonably definite direction. We can also de- 
fine energy bins that further distinguish these regions, in 
order to study energy-spectrum-flavor connections. We 
define the operators 



P, 



n v V 



id{a} 



(6) 



Then we can write the equations of motion Q as 

jP a = \ a B xP a - g P a x ^(1 - cos0 a , h )P b (7) 

a, b 

where g = y/2n v GF and x a ^ — cos6> a .b is the cosine 
of the angle between the momentum directions in the 
two groups. The angle-flavor-energy joint distributions 
of neutrino density are now determined by the P a . For 
example, the respective neutrino densities in the two fla- 
vor states r and s are given by 



n r -n, y (l + ^Pi 3 ))/2, 

a 

n s = n, y (l-5>i 3 ))/2. 



(8) 



We note that in the isotropic case of assumption (b), 
we can make the replacement — x a,b)Pb - * J where 

J = J2 a ^ a - I n this form the equations (0 have been 
studied extensively, both as a curiosity and as an essen- 
tial tool in studying the evolution of the early universe 
in the case of large initial neutrino-anti-neutrino asym- 
metries. The qualitative outcome, in this isotropic case, 
can be stated fairly simply: If we start with the individ- 
ual neutrinos in flavor eigenstates, then the timescale for 
evolution of the macroscopic properties is of order Tj,, 
just as it would be in the absence of the neutral current 
interactions. There is a much discussed effect of the neu- 
tral current interactions during this time period, namely 
the synchronization of oscillation frequencies. But this 
synchronization does not change the general evolution 
time-scale for any distribution of neutrinos that is rea- 
sonably clustered in energy and that begins from a flavor 
eigenstate. 

In the present paper we consider some solutions of Q 
that begin with non-isotropic distributions. We find cases 
that give fast evolution of macroscopic properties, in con- 
strast to the behavior described above. The non-linear 



equations © for the non-isotropic case with unequal dij 
are so complex that we have little understanding of the 
full range of possibilities. However, below we will de- 
velop a series of special cases, which perhaps can serve 
as a basis for some educated guesses about the broader 
issues. In each of these cases we seek interesting behavior 
involving macroscopic changes, but on the fast time-scale 
Ts- We shall also drop the first term on the RHS of 
in favor of an initial condition corresponding to the first 
(flavor-mixing) term having been turned on (by itself) for 
a brief period time, leading to a tiny rotation in flavor 
space, so that we no longer have quite an eigenstate of 
the original flavor operators. We explain later why this 
is an illuminating way to proceed. 



2. TWO GROUPS OF NEUTRINOS WITH 
STANDARD COUPLINGS 

By a group, we mean an assemblage of neutrinos, each 
initially in the identical flavor state and all moving in the 
same direction, in a "laboratory" frame in which we do 
the calculation. A distribution of energies within a group 
is allowed. When we have two groups, A and B, we can 
assemble the operators for the first into a collective oper- 
ator t = ®i where the notation indicates a sum only 
over group A, and a collective operator £ = J2b The 
Hamiltonian p|l. discarding Hr, and taking the contri- 
bution from the neutral current part only is now, 



Him = 



G F (l - x a ,b) 
y/2V 



(9) 



Note that the angular factors, e.g. 1 — xa.a — 0, elim- 
inated all interactions among the neutrinos within a 
group. As explained above, we will first use the equa- 
tions for the operator expectations Pa = (t) /(n v V), 
Pb = (C)/(n u V) as in Q, viz, 



^r?A = --r?B = -g(l - x a ,b)Pa x P B ■ 
at at 



We take initial conditions such that pfl — cab and 



(10) 



P 



(1),(2) 



0, P J 



(1) 



e p(2) 

ei P B 



€2, where the coefficients 



ca,b are of either sign and of order unity, and where ei,2 
are very small. Then we ask whether there are significant 
changes of occupancies Pab m the shorter time scale Tg. 

The system (fTOf) conserves -Pa(*)+-Pb(*) = a, so that (fTOf) 
can be written as, 



dt 



P A {t) = -g(l-x AlB )P A (t) x 5, 



with the initial condition p[ l \o) = c A Si t s. We see that 
P A precesses around the nearly parallel vector a with 
a rate that is of order eT^ 1 . Thus by our definition of 
short-term, the short-term change in P is small because 



(11) 
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of the factor e. In the case in which \ca — Cb\ >> e the 
change in P is doubly inconsequential, since then Pa is 
almost parallel to a and the amplitude that oscillates is 
also small. 

3. TWO GROUPS OF NEUTRINOS WITH 
SLIGHTLY ALTERED COUPLINGS 

If we slightly alter the Hamiltonian for the above case 
in a way that destroys the SU2 symmetry in the internal 
space we find an interesting generic phenomenon, one 
which will serve as the prototype for similar phenomena 
in cases with three or more groups that interact through 
the correct neutral current coupling. We take 

HioI = ^^^V ^ + t (2) C (2) + 7 r (3) C (3) ) , 

(12) 

where the SU2 symmetry is broken when 7 / 1. We 
write the equations of motion for the P's in terms of a 
rescaled time coordinate, t' = g[l — XA,B]t, 

d p(l) _ p(3)p(2) -,p(2)p(3) 

d p(2) _ p(3) p (l) , -,p(l)p(3) 
~j^r A - ~^A B + 7-TA r B , 

d p(3) _ p(2) p (l) p(X) p (2) ( ,~s 
~^j^A '-A B ~ "A r 5 ) \ lA ) 

with three more equations generated by A <-> B. Intro- 
ducing the new variables, 

x = 2Pf ; y =(Pi 1) ) 2 + (Pf) 2 + (Pi 1) ) 2 + (Pf) 2 , 

- _ o/p(2)p(l) p(l)p(2)s 
z — A \ r A B A B 

W = 2{pfpf+Pfpf), (14) 

we consider only solutions in which the conserved quan- 
tity, Pf ] + P B 3) = G. For this case O yields the follow- 
ing closed set of equations, where x = (d/dt')x, etc., 

x = 2z ; y = —2xz, 
z = — 2x(y + jw) ; w = 2 r yxz. (15) 

Once two constants of motion are set, 

x 2 

y + w/7 = ci ; y = -y +c 2 , (16) 
we obtain the second order equation, 

l = x 3 (1 _ 7 2 } _ x[2c2{l _ 7 2 } + ^2^] (17) 

We take an initial configuration in which x = 1, y = 
e with e << 1, 2 = 0, w = 0, giving c\ = e, C2 = 
1/2 + e, i(0) = 0. This corresponds, for example, to 
beginning with group B entirely in the s flavor state, 



P 




FIG. 1: The flavor turnover of group A, from the solution 
of eas. 1131 plotted against the scaled time (t' in text). P 
represents 2P^ in text, with P = 1 indicating 100% r flavor 
for group A and P = — 1 representing 100% s flavor. The 
parameter e is set at .0001, with curves plotted for values 
7 = 0, .3, .6, .9, the heavier lines standing for the lower values 
of 7. When 7 > 1.0 the curve is indistinguishable from the 
line P = 1. 

and with group A in the state derived from the r flavor 
state by rotation through an angle w \fi. We show in 
fig. 1 the results of the solution of (|17|l for the case of 
e = .0001 for various values of the parameter 7. We 
see the complete and abrupt trade of flavors between the 
A group of states and the B group of states at a time 
that is nearly independent of 7 for the range < 7 < 
.7. As 7 approaches unity the turnover point starts to 
recede. For 7 > 1 the turnover phenomena disappears, 
the flavors remaining stuck essentially on their original 
values over periods of scaled time t' < e _1 , as indicated 
by the solution (|ll|l for this case. 

Qualitatively one can understand this turnover behav- 
ior, for example in the case 7 = 0, by observing that 
when e = 0, (|17J) becomes just x/2 = x 3 — x, which has 
the familiar kink solution x = tanh(t — to) (as well as the 
solution x = ±1). Our above initial values for the case 
of very small e then tie to a periodic array of well- sep- 
arated up-kinks and down-kinks, centered at the points 
t n . Locally, each kink is nearly of the form tanh(i — t n ). 
However, we are hear interested only in determining the 
time elapsed up until the first turnover. In fig. 2 we show 
the time development over roughly this shorter span, for 
a series of values of e's equally spaced in (— loge) as the 
values of e are successively reduced. 

From the equal spacing of the plots we clearly see 
that the transition time to the first turnover increases 
as (—loge). We can understand this from the f(t) = 
tanh(< — to) solution, where if we fit a boundary condi- 
tion /(0) = 1 — e, we find t = — loge. 

Looking at (|17|l the reader might say, "Of course some- 
thing drastic will happen as we pass through the value 
7=1 from below to above. The potential in the equation 
now goes to —00 as x — > ±00. The underlying field the- 
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P 




FIG. 2: As in fig.l, now showing the set of curves for 7 = 0, 
with values e = .00001 times 1,2,4,8,16,32,64. The curves 
farther to the right correspond to the smaller values of e. The 
equal spacing is indicative of the (— log e) behavior that is one 
of the principal results of this paper. 



where the constant depends on the initial conditions. 
In fig. 3 we show numerical plots for the time evo- 

(3) 

lution of P\ when the four groups consist of two that 
are respectively in the ±z directions, and two others that 
are in opposed directions at a very small angle to the ±z 
directions. The initial assignment of flavor for the four 
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ory must lack a ground state." But we should emphasize 
that in the context in which l|17|) was derived, first, the 
input physics always limits values to numbers determined 
by the neutrino density. 

FOUR GROUPS OF NEUTRINOS WITH 
STANDARD COUPLINGS 

The last exercise was interesting, although not directly 
applicable to neutrinos themselves, since for the neutrino 
case we have 7=1, with an SU2 invariance in the fla- 
vor space. But since this case lies on the boundary be- 
tween the stable 1 and the unstable case, it is not sur- 
prising that with a more complicated initial configura- 
tion we find parameter regions in which the instability is 
present when we use the usual neutral current couplings. 
To illustrate we consider the equation for the densities, 
Pa,Pb,Pc,Pd, for the case of four groups. In the two 
group case, the angular dependence of the A — B interac- 
tion was subsumed in an effective coupling constant, g. 
In the four group case we define six effective couplings 
9a.b = V2(l — XA,B)n v Gp, etc. The equations of evolu- 
tion are then, 

d - - - 

—P A = ~9a,bPa x P B - g A ,cPA x Pc 
dt 

-9a.dPa x Pd , (18) 

plus six more equations for (d/dt)PB and (d/dt)Pc, ob- 
tained by permuting subscripts. The functions Pd can 
be eliminated through, 

Pd + Pa + Pb + Pc = const. , (19) 



i. e., stable in our defined "short term" sense over the time period 
of order Tg. 



FIG. 3: The flavor turnover of group A, from the solution 
of eqns. 1181 . where we have taken the cosines of the an- 
gles between the groups as, xa,b = —1,xa,c = .9, xa,d = 
— .9,xb,c = —-9,xb,d ~ .9,xc,d = — 1. Initial flavors 
were taken as r for groups A and C, and as s for groups 
B and D. Then a small rotation of angle e around the (2) 
axis in flavor space was performed for the initial conditions 
of group B only. Shown are the curves of P = P A for values 
e = .01, .001, .0001, .00001, .000001, .0000001. 

groups consists of two sets of flavor r and two of flavor 
s, are explained in more detail in the captions. As in the 
examples in the previous section, we took some flavor 
mixing, of order e, in the initial state, and plotted curves 
for equally spaced values of (— loge). 

We see a result very much like that of the previous 
section for the case of 7 = 0. The time required to reach 
the first turn-over, in which (P A ) reaches the negative 
of its original value appears again to be proportional to 
(—loge). Now, however, we are dealing with the correct 
coupling scheme, with its full rotational symmetry in the 
internal space. We have explored the space of input pa- 
rameters (in numerical calculations) sufficiently to say 
that there exists a substantial region of parameter space 
in which the behavior is generically like that shown in 
figs. 3, as well as a region of parameter space in which 
there is only very tiny change in occupancies over the 
time interval shown in the figures. In the latter cases 
the changes remain quadratic in e for small e . We do 
not have an analytic reduction of the equations that pro- 
vides us with the function of the parameters that defines 
the boundary of the region, as we had in the broken- 
symmetry two group case parameterized with 7. From 
the perspective of the long time-scale, that is to say, the 
usual neutrino oscillation time-scale, we would describe 
the cases where anything significant happens over the 
time t <<Tl as being unstable. 
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Looking more closely at the possible physics in the sit- 
uation that we have simulated, the reader who puzzles 
through the angular domains described in the caption to 
fig. 3, will see that we began with two groups moving in 
nearly the same direction (say the z direction) , one with 
flavor r, the other with flavor s, and two more groups 
moving nearly in the opposite direction, also with one 
flavor r and the other flavor s. Then the reader might 
ask, "If the alleged turnover just changes all r's to s's 
and vice versa, then it hasn't done much, has it? We still 
have a group of r's and a group of s's moving more or less 
in the z direction, and a group of r's and a group of s's 
moving roughly in the —z direction as well." However, if 
the energy spectra of the original r groups was greatly 
different from that of the original s groups, then we have 
indeed made a major change, since after the turnover, the 
s flavors have the energy spectrum of the original r fla- 
vors, and vice versa. We also note that by choosing other 
geometries and parameters we can contrive situations in 
which there is gross change in the angular distribution 
for a a particular flavor (though, of course, absolutely no 
change in the total angular distribution for neutrinos). 

Pending resolution of questions about the behavior 
of systems with more continuous distributions in angle, 
which we discuss a little more in the next section, we 
shall not address the complexities of application of these 
ideas to the environment of the supernova, where there is 
much anisotropy. However we should point out that the 
possible interesting unstable behavior can be expected to 
occur on a faster time scale, and farther into the super- 
nova core, than in a theory with only istropic distribu- 
tions. We note that the calculations presented in ref. [25| . 
where a flux- weighted angular average of (1 — cos 6^) is 
used, do not capture the physics of the present paper. 
Though it takes note of the angular factor, this approach 
reduces the evolution equation exactly to the form of the 
isotropic case with an effective coupling constant, and 
therefore removes the possibility of finding the rapidly 
growing modes that are the focus of the present paper. 
Moreover, it is in the vicinity of the neutrino-sphere that 
the angular distributions of the electron neutrinos and 
of the mu or tau neutrinos are very different, giving, it 
would appear, the maximum possibility for rapid spectral 
and angular exchanges to occur. 

A further caveat for the supernova case is that we have 
absolutely no insight into the nature of our instabilities 
in a world in which three flavors of neutrinos are present. 
This question calls for further examination. 

5. LINEAR STABILITY ANALYSIS WHEN THE 
INITIAL STATES ARE FLAVOR EIGENSTATES. 

We have presented two examples in which a small 
flavor-mixing perturbation in the initial condition leads 
to exponential growth followed by nonlinear oscillations 



at early times. In these examples, the neutrino distribu- 
tion in momentum space was rather artificially limited to 
a small number of discrete directions in this space. Ide- 
ally we would now examine the behaviors that are pos- 
sible when the initial flavor distribution in momentum 
space is more continuous. Direct numerical simulations 
over the full range of the non-linear oscillations arc ruled 
out by their complexity. But we note that in the discrete- 
angle examples of the last sections the early exponential 
growth period could have been studied in the approxima- 
tion of keeping only the linear response to perturbations 
of the initial conditions. We have not even been able 
to carry this programme through in any general way for 
the case of continuous distributions in angle, since the 
non-linearities defeat the partial wave expansion that is 
a natural beginning point. However we can use a partial 
wave expansion for the case in which the unperturbed 
solution is isotropic, in order to confirm our belief that 
there are in this case no growing modes of any multipo- 
larity. 

We begin with a set of neutrino density functions Q(£l) 
defined for each direction of space (i.e. the momentum 
space of the neutrinos), 

Q(n)dn = -L- V Si, (20) 

n v V 

icfdn} 

where the notation indicates that the sum is to be taken 
over all states for which the direction of momentum lies 
within the element of solid angle d£l. We introduce a sec- 
ond set of densities, R(Q), for some other neutrino states 
in the ensemble, for which the one-particle operators are 
labeled as oV , 

R(Q) dQ = — V S v . (21) 
n v V — ' 

i'Cfdfi} 

The physical distinction between Q(£l) and R(Q) could 
be that they represent different, disjoint parts of the en- 
ergy spectrum (which could have differing initial flavor 
occupancies in some applications.) The notation i' would 
then indicate that the sum was over a set of energy states 
disjoint from those indexed with i. Alternatively, the 
R(£l) variables could represent the densities of antineu- 
trinos in which case i' designates the operators for an- 
tiparticles. The way in which anti-particle densities enter 
the equations which we shall write has been treated in 
detail in ref. > an d m several other references already 
cited, although with some variations in the definitions 
of the density operators. We use the negative of the 
density matrix elements defined for antineutrinos in ref 
[lij . making the resulting structure completely symmet- 
rical between neutrino and antineutrino (see Appendix 
A). The commutation rules, both for the case in which 
Q and R are distinguished by spectrum and in the case 
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where they are distinguished by particle- antiparticle clas- 
sification, are then, 

[Qi(n) , Qjin')} = -\ e«,**(n - - ( 22 ) 



and 



[^(fi),^')] 



2i 



^ ei>J - fc 5(0-0')^(n) , (23) 



and the equations of motion, in complete analogy to (JJJ 
are now, 

^Q(Q,t) = gQ(n,t)x{ J dQ.'[Q(n',t) 

+5(fi / ,t)](l-cos0 n ,n')} , (24) 

and 

— 5(n, t) = ^(n,t)x{ /" dfi'[Q(n',t) 

+ J R(O',*)](l-cos0n,n')} • (25) 

We first consider the application of \2ty and (|25[) to a 
system that is totally isotropic (in momentum space) in 
its unperturbed state, as well as being initially in a flavor 
eigenstate. We write the two densities as 

Q(Q,t) = Q (0) +X) y i,m(M)gU(t) , 

I, in 

R(Q,t) = R (0) + J2 Yl 

,m(6,4>)n,m(t) , (26) 

I. in 

where we take space-isotropic and flavor-diagonal values 
for the unperturbed densities, 



<3(o) = ak ; i? (0) = bk 



(27) 



with k a unit vector in the (3) direction in the internal 
space, The perturbations have been expanded in spheri- 
cal harmonics (in the momentum space) with coefficients 
<fz,m,n,m- With these choices, the equations j^H an d 
(|2*5|l are obeyed if we set qi(t) — and ri{t) = 0. Thus, 
the linearized equations for the time dependence of the 
perturbations are homogeneous, with time independent 
coefficients. They follow directly from substitution of 
into J2U and J2HJl, yielding, for the cases I = 0, and 
1 = 1, 



i d , . d _^ . . 

= [-afo t o(t) + &<fo,o(*)] x k , 



(28) 



9 l i^,m{t) = [(y + b)qi, m (t) + af 1>m (t)} x k, (29) 



g l -Tjfi, m {t) = [(y + a)ri, m (t) + x * • (30) 

For values I > 1 we obtain. 



9 -n%m(t) = {a + b)qi m (t) x k 
at 

g^^in.^t) = {a + b)f hm (t) x k . 
at 



(31) 



It is easy to see that in this case, with isotropic unper- 
turbed distributions, all of the eigen-frequencies of l|29() - 
(|31|l are real. Note that the complication due to the mul- 
tiplicity the vector components, q, f in the above equa- 
tions is superficial. The eigen-frequencies of the system 
11291 , 13U|) are given exactly by the eigenvalues of the ma- 
trix, 
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b 



(32) 



If we had introduced angle dependence into the zero'th 
order distributions l|27|l . the linearized equations for the 
perturbations would still have time independent coeffi- 
cients, but all values of I would have been coupled to- 
gether. We therefore arrive at no conclusion as to the 
presence of growing modes in this case. 



6. QUANTUM CASE 

We review the development to the present point. The 
original Hamiltonian for several groups of neutrinos is 
a sum of pairwise interactions among collective opera- 
tors for each group. In the Heisenberg equation of mo- 
tion for these operators we then replaced the operators 
by their expectation values, and we replaced products of 
the operators by products of expectation values, follow- 
ing assumption (a) of sec.l. We considered only initial 
states for which this assumption is true, but the assump- 
tion cannot remain exactly true as time progresses. The 
applications of the classical equations proposed in the 
literature demand that this classical approximation be 
good over the longer timescale Tj,. A heuristic justifica- 
tion might be that when the number of particles in each 
group is large, then the collective operators become, in 
effect, the operators for a system with very large spin. 
We are accustomed to the phenomenon of large quan- 
tum numbers leading to classical behavior. One of our 
aims in the present paper was to test the classical ap- 
proach by comparing with the direct numerical solution 
of the models, the latter requiring up to several hundred 
spins, in order to estimate the dependence on the particle 
number. Here we address this question through numeri- 
cal solution of the complete Schrodinger equation, for the 
the simpler configurations considered above. In appendix 
B, using the results of this section, we shall further illu- 
minate the use of the term "quantum" in distinguishing 
this approach from that of the previous sections. 
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We begin with the simplest interesting example, the 
two group model (|12fl , where the internal SU2 symmetry 
is broken by taking 7 = 0, but where we now start with 
an exact flavor eigenstate, with the parameter of the pre- 
vious calculation, e, taken as zero. In the initial state the 
neutrinos in the upper group all carry the r flavor and 
the neutrinos in the lower group all carry the s flavor. 
Results are shown in fig. 4. 



0.5 




-0.5 



FIG. 4: Curves for the model of sec. 3, but now based on 
a complete quantum mechanical solution using the Hamil- 
tonian, and with no "seed" of flavor mixing in the initial 
condition. Now the turn-over effect comes about quantum 
mechanically and is dependent on the particle number N. 
We plot P as a function of scaled time, for values N = 
8, 16, 32, 64, 128, 256, 512. The curves are to be compared to 
those of fig. 2. 

Comparing with the plot of fig. 2, which showed the 
time development for the classical case, but with e ^ 
we find a surprising similarity. We can state the result 
of the comparison as follows: in the classical case, seeded 
with an initial perturbation that is of order e, the time for 
flavor turn-over is of order (— Tgloge); in the quantum 
case 2 , with no seeding of initial flavor mixing, the turn- 
over time is of order T$ log(iV), where N — n v V is the 
number of neutrinos. The shape of the turn-over curves 
is almost identical. 

Carrying out the same calculation for the case 7 7^ 
leads to a family of curves much like those shown in fig. 
2, again with the logiV — > (—loge) correspondence. In 
the solutions of the complete problem based on i|12f) the 
dependence on N is twofold; first the implicit depen- 
dence through the coupling parameter g = 2\f2Gn v and 
second, the residual explicit dependence that enters the 
above result. It is this residual dependence that shows 
the logartihmic behavior in TV for large values of N. 



2 As we saw, N is not an explicit parameter in the equations for 
the P's in the classical case; it was removed from the equations of 
motion when the intensive density operators P a and the coupling 
parameter g were introduced. 



We have also looked at the direct solution of the evolu- 
tion under influence of the complete forward Hamiltonian 
for the example of four neutrino groups that was treated 
classically in sec. 4. The results are completely in accord 
with our expectations derived from the above, namely, 
turn-over times that are of order T^logiV, and the rela- 
tion (— loge) <-> log N. We cannot present quite as good 
curves of the N dependence in this case because of com- 
putational constraints. The method of direct solution 
takes advantage of the fact that each group of neutrinos 
is represented by total spin operators for that group; thus 
with a given total number of neutrinos the dimensional- 
ity of the calculation goes up with the number of groups. 
Still the calculation lends support to the conclusion that 
the relation of classical instability to quantum speed-up 
is generic. 

In applications to supernova neutrinos we would usu- 
ally expect that (— loge) >> log7V e ff, where we estimate 
e « T s 5m 2 /E, and N eff is the number of neutrinos in 
a volume of dimension cTs. Thus the classical approach 
based on (fTTTf) would appear to be sufficient in the su- 
pernova application. The spontaneous effect, i.e. the 
transformations that the neutral current interactions can 
bring about in the absence of a flavor-mixing term, in the 
absence of flavor mixing in the initial state, will be small 
compared to the seeding effects of the mass 2 terms, which 
produce the small mixing of order e that in turn triggers 
growth of the classically unstable mode. It would be 
worth confirming this conclusion in detailed simulations, 
but that is beyond the scope of the present paper. 

We mention that there exists a system which, in effect, 
is governed by a model that is isomorphic to the model 
of sec. 3 with parameter 7 = 0, namely the case of the 
encounter of two photon clouds, where one or both of the 
clouds has large circular polarization, and where the in- 
teraction is the given by the standard Heisenberg-Euler 
formula. In ref. 26] it is shown that large polariza- 
tion exchange can occur over times that are many orders 
of magnitude smaller than one would estimate from or- 
dinary photon-photon cross-sections, in analogy to the 
rapid flavor exchange that we find in the above models. 

We are now in a position to state a preliminary conclu- 
sion with respect to assumption (a) discussed in sec. 1. 
We indeed have evidence that when the number of par- 
ticles approaches 00, the expectation value factorization, 
generically xaj) = (<x,) x (<Jj), if present initially, re- 
mains true for all times on the short time scale. However, 
this statement holds only up to the (log TV) -1 corrections 
to reaction rates. 

These relations between the numerical calculations 
from the complete Hamiltonian and the results in the 
"classical" approximation can be elucidated in our sim- 
plest two-group case in an analytic approach. Again 
choosing the model of sec. 3, with the Hamiltonian l|12(l 
and taking 7 = 0, we introduce a set of operators that 
are bilinear in the variables r and £, 
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x = iT (+) C (-) . U = iT (-)^+) ; J / = C ( + ) C ( ~ ) , 

z = r'-'rW ; w = r< 3 > , (33) 

where = + ir^ 2 ' etc. 

We now write Heisenberg equations of motion for these 
operators by taking commutators with Hf OI of l|12(l , spe- 
cializing to the case, xa,b = 0, and taking the con- 
served quantity + £( 3 ) = 0, We measure time in units 
Ts = (Gfii„) _1 and obtain the closed set, 

Nx = w(z + y) + w 2 ; u = —x ; Ny — wx — uw, 

Nz = xw — wu ; Nw = —x + u . (34) 

Now we take (|34l) as c-number equations for the expecta- 
tion values of x, y, z,u,w . We choose the initial system 
to have N states in group A all occupied with neutri- 
nos of flavor r, and N states in group B, all occupied by 
neutrinos of flavor s. Then the initial values of the vari- 
ables defined above are x = y = z = u = and w = N. 
For this case the set l|34l) leads to a single second-order 
equation for w = w/N, 

r] 2 2in 2 

|j<D = 2w(w* - 1) - fj?L . (35) 

where the initial condition is now w — 1. In appendix 
B we give a fuller explanation of the above steps, and in 
addition, we restore U to the calculation to show the sense 
in which the last term on the RHS in (|35J) is a quantum 
correction to the classical approximation. 

If, instead of using the equations l(3*4")l for the bilinear 
forms, we had used the Heisenberg equations for f and 
for C taken now as classical quantities, and only then had 
introduced the variables x, y, z, u, w as defined above, we 
would have retrieved the set l|34ll with one modification; 
the first equation would now be just, Nx — w(z + y), 
and would be replaced by -j^w = 2w(w 2 — I). This 
set now lacks the term that drives the flavor turnover in 
time Ts \og(N). Indeed, as remarked earlier, it has no N 
dependence. 

We can derive the logarithmic behavior analytically 
from l|35|l . capitalizing on the fact that when N — > oo, 
the solution is the familiar w = tanh[(t— to)] kink solution 
in a A0 4 theory in one dimension, then showing that for 
large N, in the time region in question, the w 2 /N term in 
(|35|l can be dropped in favor of changing the initial value 
of w to I — 2/N, this in turn determining to = log(2iV). 

7. DISCUSSION 

We began with the standard classical equations of 
motion for neutrino flavor densities under conditions in 
which the time-scale of the neutral current interaction Ts 
is short compared to the flavor oscillation time-scale Tf,. 



Our first object was to examine solutions of these clas- 
sical equations for the case of initial conditions that are 
flavor-diagonal, or nearly so, but with flavor-dependent 
spatial anisotropies of the initial momentum distribu- 
tions of the neutrinos. We then solved for the time 
scales for significant change in the distributions. Al- 
though our particular solutions are periodic, in any phys- 
ical context with randomized initial conditions these time 
scales would be indicative of the time required for flavor- 
equilibration. 

In our actual calculations based on the classical equa- 
tions we discarded the neutrino mass term (which gen- 
erates conventional oscillations) in favor of initial condi- 
tions that depend on a parameter e which is to be inter- 
preted as the amount of neutrino mixing that develops, 
starting from a flavor-diagonal state, on the short time- 
scale Ts- 

Depending on the details of the initial distributions, we 
found two dramatically different possible behaviors, sep- 
arated by paper-thin transition zones. In one case there 
is essentially no evolution of occupancies over the short 
time scale Ts- In the second case there is large evolution 
on the time-scale (— Is log e). The simplest model that 
displays this transition is the two-group model of sec. 3 
where we changed the neutral current coupling somewhat 
by introducing the parameter 7, and where the physical 
value 7=1 turns out to be exactly the boundary of 
short-term stability. But in the multigroup cases that we 
studied, there are regions of instability for the physical 
case 7 = 1, as well. 

Our second objective in this paper was to better under- 
stand the limitations of the classical approximation that 
underlies almost all of the literature on this subject. We 
recall that our path to these equations involved two steps: 
I) Replacement of the complete neutral-current Hamilto- 
nian by a forward Hamiltonian that includes only the mo- 
mentum modes occupied (by neutrinos of either flavor) 
in the initial state ; 2) The factorization assumption JSJ. 
In sec. 6 we investigated this latter assumption by direct 
computations of evolution beginning from the forward 
Hamiltonian. In order to achieve the maximum clarity 
in the interpretation of results, we began from exact fla- 
vor eigenstates, e = 0, in carrying out these solutions. 
For this case, of course, the classical equations give no 
evolution. Under conditions for which the classical equa- 
tions are stable, in the sense defined above, the numerical 
solutions for the quantum case, indeed show only a tiny 
evolution over the time scale Ts, by an amount that de- 
creases (as iV _1 ) as is increased. 

However when we look at the same comparison under 
conditions in which the classical solutions are unstable, 
the results are very different. In this case the numerical 
solutions for finite iV can undergo complete turnover of 
flavors, from one part of momentum space to another, in 
a time that is apparently of order Ts log N. The inter- 
esting comparison is now to the classical approximation, 
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now with a small value of e, where we found similar devel- 
opment over the time scale (— Tgloge). These observa- 
tions are based on simulations that are, for all cases but 
one, a bit inadequate in covering the range in N required 
to make a firm conclusion as to the logarithmic behavior; 
hence our use of the word, "apparently" , in the above. 

In the simplest model that we considered, that of sec- 
tion 3, we were able both to calculate for large enough 
values of N to make the above conclusion fairly firm, and 
also to find an analytic approach that supports the log- 
arithmic behavior. Qualitatively we could describe this 
latter approach as being the next step in a heirarchy of 
approximations that begins with the classical approxi- 
mation. The classical approach writes equations for den- 
sities, quantities bilinear in the neutrino fields, with a 
factorization assumption to close the set. This is now 
replaced with equations for quantities that arc bilinear 
in the densities, again with a factorization assumption to 
close the set. The solutions of this second set of equa- 
tions well replicate the results of the numerical study of 
the quantum case over our whole range of parameters. 

From these results, and further simulations that we 
have not presented here, we have what we believe is a 
generic relationship between the instabilities of the clas- 
sical approximation and the T$ log N behavior of the full 
models. Then returning to item # 4 on the list of pos- 
sible neutrino collective phenomena at the beginning of 
this paper, we have at least a partial answer to the ques- 
tion of whether there is ever a "speed-up" of the sort 
proposed in ref. 0] but questioned in ref. 17]. Ref. 
|l5| addressed the case with an initial condition in which 
the two neutrino flavors occupied different spectral re- 
gions but in which the angles between the momenta were 
distributed over the whole region — 2n, and asked for 
the characteristic time for the spectral rearrangement. 
The conjecture presented was that the evolution would 
be at the "speeded-up" rate characterized by T5. Be- 
cause of the scatter in couplings, as well as inadequate 
numerical methods, ref. [lfi| achieved simulations with a 
maximum value of only N— 14, far too small a value to 
find the logarithm if that were the answer. In the present 
paper, we have been able to study cases with much higher 
values of N through taking a few beams narrow in solid 
angle, and we find the speed-up discussed above, for cer- 
tain initial conditions. 3 Ref. ^3] took all couplings to 
be the same, as for the totally isotropic "classical" equa- 
tion, and solved the full problem analytically, finding no 
"speed-up" . Since the classical isotropic case is stable, as 
we saw in detail in sec. 5, the lack of a speed-up in the 
full solution to the model of ref. ^3] model presents no 



3 Note that even with the log N factor included we can legitimately 
claim a "speed-up" in many circumstances, since T%, is often 
many orders of magnitude greater than Tg . 



inconsistency with our generic relation between stability 
of the classical case and speeded-up behavior for the full 
model. 



There is a difficult question that we have not addressed 
here. We have studied systems that consist of plane 
waves that fill our quantization volume. In our mod- 
els the participating particles lie on top of each other 
for the whole time interval over which we study the sys- 
tem. In actual physical systems that we might consider, 
the particles moving in a given direction must have some 
characteristic length that determines how long they stay 
in contact with those moving in the opposite direction. 
We are presently not in a position to say exactly how this 
affects the analysis of the "quantum" evolution effects in 
any realistic situation (we believe, as do other authors, 
that it affects the "classical" evolution not at all). We 
can try to gain insight with calculations which are sim- 
ilar to those reported in sec. 6, but where instead of 
having two groups of particles in contact with each other 
over our complete time interval we focus on the state 
of a single set, moving in the same direction, and take 
this set to be serially in contact with a number of groups 
of other particles, one group at a time. We can state 
some prelimary results of these calculations in a negative 
form, first by posing the question, "As our test group is 
affected first by one other group, and then by another, 
and by another, is it just a question of determining, af- 
ter each encounter, the probabilities of transitions in the 
test packet, and then proceding to the next encounter 
with an flavor-diagonal ensemble embodying these prob- 
abilities?" If the answer were yes, then of course there 
would be no evolution beyond that which is be achieved 
in a single short encounter time, multiplied by the num- 
ber of interactions. Based on simulations that we have 
done, reinforced with some analytical considerations, we 
can state firmly that the answer to the above question 
is "no". There is a coherent effect, in spite of the inde- 
pendence of the colliding packets, which gives rate en- 
hancements. We shall return to this subject in a future 
publication. 



There are also possible questions relating to our use 
of the "forward" Hamiltonian (JJJ rather than the full 
Hamiltonian for the system; that is, to the neglect of all 
but forward scattering processes. Here we believe that 
we are on a firm footing, as long as we consider times 
that are very much less that T scat . Inclusion of the non- 
forward parts gives a forest of oscillating terms that adds 
up essentially to nothing over this time scale. We know 
of no reference that treats such questions in detail, but 
we do note that the whole literature on what we have 
called the classical approach is built on this assumption. 
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APPENDIX A. INCLUSION OF 
ANTIPARTICLES. 

Since the above work did not explicitly deal with 
neutrino-anti-neutrino mixtures, and such mixtures will 
enter any application to which the formalism is relevant, 
we show how anti-particles fit into the general framework 
without changing the structure of the equations. We de- 
fine the operators, 

a\ +) = {a^a* + 

ctH = (ajjtaj + (ajjtoj, (36) 

where now the barred operators create and annihilate 
the anti-neutrino states for flavor (r, s) and momentum 
Pj. We note the relation, 

= °i 3) = K) f s r - ( a PS s 

-(aj)ta5 + (a|)t55 . (37) 

Now forming the combinations, erj 1 '' = (<xj + <r^ 

<rj 2 ' = — i(crj + ^ — <r j ^) for each mode, and substituting 
in Q , where the index j labels momentum only, we ob- 
tain the complete relevant part of the (neutral-current 
plus oscillation) forward Hamiltonian, including antipar- 
ticle operators, as can be shown by a tedious calculation. 
By relevant we mean that all the remaining interaction 
terms, as well as the kinetic energy terms, commute with 
each of the 5\. The <7; operators, as extended above, still 
obey angular momentum commutation rules, so that all 
of our previous formalism for following the distributions 
in time is applicable to groups that contain mixtures of 
particles and antiparticles, with the provision, of course, 
that the total amount of one particular flavor is counted 
by subtracting the number of anti-particles of that flavor 
from the number of particles, as in Ij37|l . Thus the equa- 
tions (0J still describe the time evolution of the system. 

We may wish to separate particles from antiparticles, 
however, putting them in disjoint groups with collective 
operators P a , where various values of the subscript a now 
both distinguish energy bins and distinguish particles 
from antiparticles, and where an anti-particle operator 
is defined, for example, as, 

Then it is easy to see that dividing the collective oper- 
ators into particle and antiparticle parts indexed with 
different a preserves the commutation rules of the under- 
lying operators, and that the equation of motions for the 
P a 's, after we make the classical factorization assumption 
(0), are again given by 10, 

-P a = \ a B xP a - gP a x ^(1 - cos9 a , b )P b . (39) 



We contrast this equation with the equations that have 
been given in |l9| and other sources, here specialized to 
the isotropic case where (in effect) we have cos0 a (, = 0, 

d„ Am 2 „ V2G F ^ _ „ 

-2p- B x -v^ x 5> - po, 

j % 

d= Am 2 V2G F = \ ^ , - „ , 

^=-^BxP J + — Pjx^-P.,) 

J i 

(40) 

where Pj stands for a particle density (or density matrix 
element) and Pj stands for an antiparticle density. 

We can reconcile (13911 with the set l|4U[l first by taking 
into account the difference in normalization of the P Q 's in 
the two cases, and the corresponding relation of coupling 
constants. 4 and then following with the replacement 
Pi — > —Pj in the equations |@DJ, from ref 0. The lat- 
ter change reflects a difference in definitions of the den- 
sities 5 Our definitions are better suited to our purposes 
because when we consider the operators of which the P's 
are expectation values, it is with our choice that both the 
particle and the antiparticle sets of underlying operators 
have the same commutation rules (angular momentum 
commutation rules modified by the factor {n^V)^ 1 ). 

APPENDIX B. ILLUSTRATION OF THE 
COMPARISON OF QUANTUM AND CLASSICAL 
CALCULATIONS. 

When we introduce h into the problem and define 
the commutation rules for the Fourier components of 
the fields in the canonical way, then forming the bilin- 
ear forms t and £ as in (0, the commutators become 
[tP,t\ 2) ] = 2i7Lrf ) ,etc. Taking X A , B = 0, 7 = 0, the 
Hamiltonian 1121) is 

tf for = ^(T«C«+T< 2 )C«). (41) 

As in sec. 6 we define 

x = lT (+) C (-) ■ „ = i T (-)£(+) ; y =C (+) C (") ) 

z = r (-V (+) ; w = t ( V t (42) 

then calculating iU~ [x,H], and obtaining (using un- 
sealed time), 

i = ^[ T (3) C (+) C (-)_ C {3) T (+) T (-)] 



4 Our normalization is such that ^ P„ ' = (n r — n s ) /n v ; that of 
ref. HI has the additional factor n u on the RHS. 

5 As noted in ref. 111! , much of the literature uses definitions that 
give equations that are even less symmetrical between par ticle 
and antiparticle than 1401 . We refer the reader to ref. [l9] for 
details. 
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^[r (3) (C (+) C H +r(-M+)) + Mr' 3 ') 2 ] 

G 

= — ^[2wz + hw 2 ] . (43) 

To get the second line we have here used + £' 3 ) = 
and evaluated one more commutator in order to get the 
product of r's in the order that agrees with the defini- 
tion of z, above. Note that some guile has been used 
to get a closed set of equations for our five quantities, 
and that the procedure is somewhat arbitrary as well, 
since we have used the noncommutative properties of 
our operators up to the point at which we can reex- 
press results in terms of our original variables. Then 
we replace them by c-numbers that we take to be their 
expectation in the medium. (Note also that the "classi- 
cal" approximation that we discussed earlier is subject 
to exactly the same criticism.) With our initial con- 
ditions x = y = u = z = and using the equations 
x = —ii, y = z (where in the latter the expectation value 
ansatz avoids an order-of-operator difference) , we obtain 
u(t) = —x{t) and y{t) = z(t). Then the two remaining 
equations of motion are easily found to be, 

vn 2 L J 

z = ^[2xw\. (44) 

We note that the last term on the RHS of the equation 
for x has the only factor of H in the three equations of 
motion, aside from the factor of h~ 2 that multiplies Gf, 
which will enter into any time-scale depending on Gf- 
Thus, in a well defined sense, this term is a quantum 
correction. Indeed, if we had calculated the equations of 
motion for these bilinear forms using the original classical 
equations for the linear forms, (|13|l for the present case, 
the result would have been exactly the above set, except 
lacking the last term in the x equation, with its explicit 
Ti. Therefore we have a basis for calling the classical 
approximation "classical" and describing the computer 
generated solutions plotted in section 6, or alternatively 
the solution of l|35|l . as "quantum mechanical. " 

For the case of our initial conditions, the reduction of 
the three equations l|43l) and l|44l) to the single equation 
(|35|l for the function w proceeds by noting first that by 
dividing the first of eqs. (|44|) by the second, and doing 
one integral we obtain, 

w 2 /2 = —z + const. . (45) 

Then, reverting to units with h = 1, scaling the time, 
choosing w(0) — TV, and defining w = w/N we obtain 
(EiJ). 
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